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We present grazing incidence reflectivity measurements in the far infrared region at tempera- 
tures above and below Tc for a series of thallium (Tl2Ba2Cu06, Tl2Ba2CaCu2 08) and bismuth 
(Bi2Sr2Cu06, Bi2Sr2CaCu2 08, and Bi2-xPbxSr2CaCu208) based cuprate superconductors. From 
the spectra, which are dominated by the c-axis phonons, longitudinal frequencies (LO) are directly 
obtained. The reflectivity curves are well fltted by a series of Lorentz oscillators. In this way the 
transverse (TO) phonon frequencies were accurately determined. On the basis of the comparative 
study of the Bi and Tl based cuprates with different number of Cu02 layers per unit cell, we suggest 
modifications of the assignment of the main oxygen modes. We compare the LO frequencies in 
Bi2Sr2CaCu208 and Tl2Ba2Ca2Cu30io obtained from intrinsic Josephson junction characteristics 
with our measurements, and explain the discrepancy in LO frequencies obtained by the two different 
methods. 

PACS numbers; 74.72.-h, 74.25. Gz, 74.25.Kc, 74.76.-w, 74.50. +r, 74.80. Dm, 78.30.-j 



I. INTRODUCTION 

The phonon spectra of high temperature supercon- 
ductors were widely studied due to the possible impor- 
tant role in the mechanism of superconductivity. Be- 
sides that, the phonons being on the same energy scale 
as the superconducting gap strongly affected many su- 
perconductivity related quantities, such as electronic re- 
sponse and the superconducting energy gap, pseudo- 
gap, Josephson plasma oscillations, and tunneling spec- 
tra. The first optical nieasurements were performed on 
polycrystalline samples .EJ The assignment of the phonon 
modes and the lattice dynamical calculations used for 
the assignment were strongly influenced by the mixed 
in-plane and out-of-plane response in these earlier data. 
As large single crystals became available, the infrared 
(IR) properties and, in particular, the c-axis phonons 
of the yttrium and lanthanum based compounds were 
thoroughly investigated, providing phonon eigenvector 



patternsBij different from the early assignment. The c- 
axis response of the bismuth and thallium classes of su- 
perconductors is known to a much lesser extent. The 
measurements done an-, these materials probed primar- 
ily th&|ft^propertiesj3~Ll or polycrystalline samples were 
usedlHJ^EJ The unambiguous information on the IR c- 
axis phonon spectra wecej-abtained on Bi2Sr2CaCu208 
(Bi2212) singlk crystalsBB and on TlsBasCaaCuaOio 
(T12223) film.EZl Extensive studies, along with the efforts 
to assign the main oxygen modes in these mnmounds 
were also performed using Raman scattering.Ej~Ej 

In this paper we present a systematic study of the 
c-axis IR phonons in Bi and Tl cuprates by means of 
grazing angle reflectivity measurements. We discuss the 
technique used and show its sensitivity to the c-axis prop- 
erties in anisotropic materials. The compaMon with the 
available data on Bi2212 by Tajima et a/Jl3 proves that 
both transverse (TO) and longitudinal (LO) c-axis fre- 
quencies can be determined correctly by this method. We 



report on the first measurements of the c-axis vibrational 
spectra in single layer Bi and Tl compound, done on sin- 
gle crystals. By comparing bismuth and thallium com- 
pounds with the different number of Cu02 layers per unit 
cell (single and double layer systems), we propose a con- 
sistent assignment of the main oxygen modes, diffeccBt, 
from those done by lattice dynamical calculations B^'Eil 
In the bismuth cuprate, using lead-doped samples, we 
identify the superstructure induced modes. The par- 
ticular sensitivity of our technique to the longitudinal 
modes enables a direct comparison of our data to the 
fine structure obseriffid_in the intrinsic Josephson junc- 
tion characteristicSjEjH which was ascribed to the LO 
phonon excitations .E^ 



II. EXPERIMENTAL 

The reflectivity measurements were performed at a 
grazing angle of 80° relative to the surface normal, with 
p-polarized light (electric vector in the plane of inci- 
dence), in the far-infrared (FIR) region (30-700 cm~^). 
Reference spectra were obtained by evaporating a gold 
film in situ on the sample surface, and absolute reflectiv- 
ities were determined by reflectivity ratio. A wire grid 
polarizer was used to select the proper polarization. The 
measurements were made using a Fourier transform spec- 
trometer (Bruker IFS 113v). A series of samples was 
measured: 

BiiSriCuOg (Bi2201), BiaSraCaCuaOg (Bi2212) and 
lead doped Bi2-xPbxSr2CaCu208 (Bi2212+Pb) single 
crystals were grown at the Lebedev Physical Institute 



by the method of re-criista 
dissolved in melted KCl.ca"E 



ization from the precursor 
As-grown samples were ob- 
tained as free-standing thin crystalline platelets and did 
not require any further treatment. The transition tem- 
peratures (and the transition widths) were Tc= 6-8 K 
(ATc = 2 K), Tc = 80 K (ATc = 2 K), and Tc = 75 K 
(ATc = 3 K), respectively. The lead contents in doped 
samples was x = 0.4, as determined by the x-ray electron 
probe microanalysis. 

Bi2Sr2CaCu2 08 single crystals of a larger size (3x4 
mm) were grown at the University of Ajasterdam by the 
travelling solvent floating zone method.E3 The transition 
temperature and the transition width were Tc = 89 K 
and ATc = 3 K, respectively. An optically clean surface 
was obtained by peeling the samples along the ab-plane. 

Tl2Ba2Cu06 (T12201) single crystals were produced at 
the Institute of Solid State Physics. The samples, with 
Tc = 82 K Epji^Tc = 13 K, were obtained by the self- 
flux methodMN 

Tl2Ba2Cu06 films came from SUNY at Buffalo. They 
were made by radio- frequency magnetcon sputtering fol- 
lowed by post-deposition annealingo The films had 
Tc = 80 K. More details on the preparation and proper- 
ties of the films can be found in Refs. ^ and |4[ 

Tl2Ba2CaCu208 (T12212) films were provided by Su- 



perconductor Technologies Inc. and-were produced using 
their standard production methodscj and had Tc = 98 K. 
More details on the electrical properties of these films can 
be found in ^. 

The samples were mounted in a conventional liquid 
He fiow optical cryostat, which allowed any temperature, 
between room temperature and 8 K, to be reached. 



III. TECHNIQUE 

The geometry of the grazing angle technique is shown 
in the inset of Fig. |l|. From the figure it is easy to see 
that s-polarized light probes the a6-plane properties of 
a crystal, whereas p-polarization gives much information 
about the c-axis response. This is a convenient method 
for extracting both the a6-plane {Sab — ^'ab^^^ab ~ ''^'ab)^ 
and the c-axis (ec = e^ + ie") components of the di- 
electric tensor, in case normal incidence measurements 
on the ac-plane are not possible due to the small size 
of samples along the c-axis, or to the mica-like nature 
of the Bi compounds. All cuprate superconductors are 
highly anisotropic in the ac, or 6c-planes. The a6-plane 
anisotropy is much smaller. It does not affect our analysis 
and is neglected in the following discussion. In this case 
the reflectivity Rp is described by the Fresnel formula for 
a uniaxial crystal: 
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where is the angle of incidence, in this case equal to 
80°. In order to fit the data, we need to assume a form 
for Sab and Sc- The most common procedure is to use the 
Drude or two-fluid model (in case of a superconductor) : 
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where the first term represents high energy interband 
transitions, the second and the third terms are supercon- 
ducting and normal carrier responses, respectively, and 
the last one is a sum of phonon contributions. In the 
normal state, the electronic background largely affects 
Sab, since the a6-planes are metallic, whereas the c-axis 
properties are insulator-like. This means that the most 
prominent phonon modes are polarized along the c-axis, 
with the afo-planc phonons being strongly screened by 
the free carriers. The reflectivity function Rp has sharp 
minima when e'^ goes to zero, i.e. at the longitudinal fre- 
quencies, while the general trend depends on £„&■ The 
phonons associated with the a5-plane are much broader 
and almost absent in the spectra, so we can assume that 
Eab is a smooth function of frequency. 



The advantage of the grazing angle technique is that 
the c-axis longitudinal frequencies are determined di- 
rectly from the measured data, without any modelling. 
From the reflectivity we can calculate a pseudo-loss func- 
tion, L{uj), defined ascZl: 
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where Uat is the in-plane complex refraction index and 
(j) = 7r/2 — arg (riab) is a weakly frequency dependent 
phase shift. L(uj) is approximately the same as the c- 
axis loss-function Im(— I/ec) (at least the maxima of 
both functions are at the same positions), and it gives us 
the frequencies of the c-axis longitudinal phonons. The 
reason for observing the LO phonon (and not TQQ fre- 
quencies in Rp is common for the c-axis plasmorO and 
phonons, and is the following: The incident electromag- 
netic radiation excites a wave inside the material, which 

has an electric field component E = Ee''^'"'~'^*\ where 
each component of E and k can be complex. If the x- 
axis lies at the intersection of the plane of incidence with 
the surface, and the z-axis is normal to the surface, then 
the x-component of the wave vector, kx, is determined 
by the incident radiation, while the z-component of the 
wave vector can be shown to depend both on the in-plane 
and out-of-plane dielectric functions, as follows: 
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where w is the frequency of incident radiation and c is 
the speed of light. Inside the medium the wave remains 
polarized in the xz-plane, and the components of electric 
field relate to each other as: 
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This illustrates the well-known fact that the electromag- 
netic wave in an anisotropic medium does not have a 
pure transverse character, but it has also a longitudinal 
component, unless it propagates along the principal axes. 
In most high-Tc compounds the c-axis plasma fre- 
quency is far below the optical phonon frequencies and 
the phonon modes are well separated in frequency, while 
the in-plane response is mostly dominated by plasma os- 
cillations of carriers in the afe-plane. Thus, the dielec- 
tric function £ah(aj) at low frequencies can be written 
as a pure plasma response, and eJyiS) near a particular 
phonon can be simplified to a single Lorentz oscillator: 
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where ujpab is the in-plane plasma frequency, lot is the 
TO phonon frequency, S is the oscillator strength, and 
Esc includes also contribution to Soo from other phonons. 
In both functions, we neglect the damping for simplicity. 



The substitution of Eqs. ^ into Eqs. ^ gives a pure 
imaginary solution for kz{uj) everywhere, except for a 
narrow frequency range, where w^ varies between uf^ — 
uj^ + Slo^/Esc and lo'^ + Suj'^/{esc — sin^ 0). In this re- 
gion kz{uj) is real. This corresponds to the propaga- 
tion of a travelling electromagnetic wave in the medium, 
which reduces the refiectivity Outside this narrow fre- 
quency range the electromagnetic wave decays exponen- 
tially near the surface, and the total reflection takes 
place. The finite scattering rate in the planes gives rise to 
the Hagen- Rubens behavior of the reflectivity, while the 
phonon damping together with the c-axis conductivity 
determine the width of the phonon lines. 



IV. RESULTS 

The room temperature spectra of Bi2201, Bi2212, 
Bi2212-HPb, T12201, and T12212 are shown in Fig. 0. 
The solid lines represent the fit, and the symbols are the 
experimental data. These results were reproduced on a 
number of Bi and Tl samples. The data were collected 
at temperatures between 8 K and 300 K. The main tem- 
perature variations in Rp occur due to the changes in 
the in-plane response. Bi compounds show more opti- 
cal phonons due to the lower symmetry. For all five 
compounds, the overall reflectivity scales with the in- 
plane response. In general, the double layer compounds 
have higher conductivity in the plane than the corre- 
sponding single layer compounds and, consequently, they 
have higher reflectivity. The optical conductivity of Tl 
cuprates usually is at least twice as high as that of Bi. 
The reflectivity Rp was fitted with Eq. |l] and Eq. ||. The 
fitting parameters: TO frequencies, oscillator strengths 
Sj, and the damping parameters 7^ of the phonons, as 
well as LO frequencies (obtained as the positions of max- 
ima in the loss- function) are listed in Table B and Table O. 
The temperature dependence of the phonon parameters 
can be seen from the comparison to the 8 K data, shown 
in brackets. The c-axis conductivities and loss-functions 
calculated from the fit are shown in Fig. g] and Fig. ||, 
respectively. It should be noted that the LO frequencies 
obtained in our experiment are not sensitive to the error 
in the absolute reflectivity and to the assumed in-plane 
conductivity, which is also fitted, while the TO frequen- 
cies depend on them both. However, this dependence is 
rather weak. In Table || we showed our values of TO 
and LO frequencies in B|i2212, together with the values 
obtained by Tajima et aZ.ll3 from the normal incidence re- 
flectivity done on the ac-plane at 6 K. In the latter case, 
the TO frequencies can be evaluated directly from the 
reflectivity spectrum. The LO frequencies are calculated 
via Kramers-Kronig transformation. From this compari- 
son, one can see that the discrepancy between the LO fre- 
quency values obtained with the two different methods, 
on different samples, is only few wavenumbers except for 
the highest frequency mode, and that the TO values dif- 



fer less than 10 cm~^. This means that ahhough we have 
a rather compHcated situation, with the mixed a6-plane 
and c-axis response, we can get quite accurate values of 
the phonon parameters: transversal and longitudinal fre- 
quencies, and consequently, oscillator strength. 



V. DISCUSSION 

A body-centered-tetragonal structure lA/mmm {D\J^) 
has been most frequently used in the interpretation of 
the infrared and Raman spectra of both bismuth and 
thallium cuprates. The unit cell for Tl2Ba2Cu06 and 
Tl2Ba2CaCu208 is shown in Fig. ^. The structure of 
Bi2Sr2Cu06 and Bi2Sr2CaCu208 resembles the structure 
of the corresponding thallium compounds, where Bi and 
Sr substitute Tl and Ba, respectively. However, the crys- 
tal lattice of superconducting Bi cuprates has orthorhom- 
bic distortions, which lead to the doubled \/2 x \/2 unit 
cell in the a6-plane. Moreover, there is a monoclinic in- 
commensurate modulation of the structure with a pe- 
riod 4.5-5 b. The factor group analysis for the space 
group li/mmm predicts 5 ^2u out-of-plane and 6 Eu in- 
plane IR active modes for single layer and 6A2U + 7Eu for 
double layer compounds, respectively. However, seven c- 
axis modes can be resolved in case of Bi2201 and Bi2212 
and only four modes are clearly seen in case of T12201 
and T12212 (Fig. |l|). Extra phonon modes are present 
in the Bi spectra due to the lower symmetry. The six 
strongest c-axis modes are observed in the Bi2212 spec- 
trum at 97, 174, 220, 328, 383, and 628 cm^^. The last, 
628 cm~^, mode is split into main and satellite struc- 
tures. The satellite at 655 cm~^ and the two additional 
phonons at 469 and 472 cm~^, which have very small 
optical strength and can be resolved only at low temper- 
atures seem to result from orthorhombic and monoclinic 
distortions. This implies that deviations from the basic 
lA/mmm structure have some, but rather small effect on 
the phonon spectrum. 

The lattice dynamical calculations for T12201, T122I2 
and Bi2212 were perfoipied, using the shell-modeSEiJ 
and the rigid-ion model.CjEil The space group M/mmm 
was used in both cases. The assignment of the modes 
and the frequencies from the lattice dynamical calcula- 
tion are given in Table III. TO frequencies were calcu- 
lated in the rigid-ion model, and TO-LO splitting from 
the shell-model. Although the shell-model is a more pow- 
erful technique, it turns out that the eigenfrequencies and 
the character of the vibrations calculated in the rigid-ion 
model are often closer to the experimental data. How- 
ever, the agreement is still not perfect. In our discussion, 
we will compare the the phonon frequencies with the cal- 
culated ones, but we will go beyond the lA/mmm group 
when we discuss the modulation induced modes. From 
now on we will indicate the oxygen in the CuO planes as 
0(1), in the SrO (BaO) planes as 0(2), and in the BiO 
(TIO) planes as 0(3), as shown in Fig. 0. 



A. Phonon modes 

Let us consider first the high frequency modes. They 
mainly involve motion of light oxygen ions. For simplicity 
we will discuss here the LO frequencies of the double 
layer compounds. The corresponding frequencies for the 
single layer compounds are given in Tables || and y, and 
the assignment of the oxygen modes is identical for both 
systems. The oxygen vibration modes are at 615cm~^ 
(628/655cm-i), 366cm~i (383cm-i) and at 328 cm^^ in 
the Bi compounds. 

We assigned the highest frequency mode, i.e. at 615 
cm~^ (628/655 cm~^), to the out-of-phase vibration of 
0(2) and 0(3), in agreement with the lattice dynamical 
calculations, given in Table IH. The splitting of this peak 
in Bi2201 and Bi2212, observable in Fig. |a,|4j3 due to 
the characteristic incommensurate modulationES of these 
compounds. This monoclinic superstructure is very often 
attributed to the incorporation of extra oxygen atoms in 
the BiO planes every 4.5-5 periods, although 0(2) va- 
cancies or element substitutions are also considered to 
explain the superstructure. The modulation period 4.5 b 
on our Bi2212 sample was measured by LEEDE3 What- 
ever the reason of the superstructure is, we use the fact 
that the modulation is not sinusoidal, but rather forms 
narrow stripes aligned along the a-axis, as indicated by 
high-resolution electron microscopyEll and x-ray detec- 
tion of higher order superstructure reflections. In this 
case, we can consider an oscillator model with every fifth 
oscillator slightly different as it could be the case in pres- 
ence of extra oxygen. Because there are four undisturbed 
and one disturbed cells in each supercell, the ratio of 
optical strength of these two modes should be around 
4, which is close to the ratio of 4.6 observed in Bi2212 
(Table ^. In Bi2201, the interaction within BiO lay- 
ers is relatively stronger, due to the weaker perovskite- 
like Cu02 block. The higher interaction strength causes 
the redistribution of the spectral weight in favor of the 
stronger mode. The disappearance of the double struc- 
ture in Pb doped Bi2212 (see Fig. IT]) supports this mode 
assignment. The substitution of 20% Bi^+ by Pb^+ re- 
laxes the strong monoclinic superstructure distortion and 
induces a weaker orthorhombic modulation with a period 
at least twice as long.E3 No other mode shows such notice- 
able splitting. Taking into account that the most signifi- 
cant displacement of ions with respect to the tetragonal 
U/mmm structure occurs in the BiO planes and, thus, 
the modulation should mostly affect interatomic bonding 
in these layers, we gave the earlier mentioned assignment 
to this mode. 

We ascribe the weak shallow structures between 500 
and 600 cm~^ to in-plane phonons. The minimum at 508 
cm^^ is also present in the a6-plane reflectivity 00312201, 
and it was reported as an in-plane phonon modeJj A sim- 
ilar observation seems to hold for the 560 cm^^ mode 
in TI22OIJ3 Small sharp minima around 470 cm"^ in 
Bi compounds are most probably zone-boundary modes. 



which become IR active due to the lattice distortion. 

The next oxygen mode at 366 cm~^ in the Tl com- 
pound (383 cm^^ in Bi) is one of the strongest modes. 
It is the most noticeable feature of the spectra in Fig. |l| 
that this mode shifts up in frequency by 61 cm^^ in Tl 
and by 44 cm^^ in Bi cuprates, going from the double to 
the single layer compound. This shift indicates that this 
mode involves the piotion of oxygen in Cu02. In the shell 
model calculationsE2lEil (Table [11 ) , the next predicted LO 
mode at 462 cm"! in T12212, and at 497 cm'^ in Bi2212, 
is too weak and too high in frequency to be assigned to 
this vibration. The calculated mode eigenvector involves 
the motion of 0(1) against 0(2) and Ca. However, the 
compound without Ca should show a lower resonant fre- 
quency, as reported in Table III. The opposite behavior 
of this mode indicates that maybe the influence of Ca 
on this mode is not so strong. If we consider the next 



available candidates for the assignment in Table [II, i.e. 
294 cm^i in T12201 and 355 cm^^ in T12212 modes, we 
again obtain the opposite frequency change from single 
layer 2201 to double layer 2212 compound. Thus, no 
mode predicted in lattice dynamical calculations could 
be ascribed to this phonon. Since the observed oscillator 
strength of this phonon is very high we think that this 
mode involves mostly the in-phase motion of 0(1) and 
0(2) against Cu. 

The mode at 328 cm~^ is clearly seen only in the Bi 
compounds and it is present in both Bi2201 and Bi2212 
at roughly the same frequency. In T12201, this mode also 
occurs at approximately the same frequency, but the opti- 
cal strength is much weaker, and the line can be resolved 
only at a very low temperature. The 328 cm~^ mode in 
T12201 shows narrowing with temperature, which is char- 
acteristic for c-axis phonon lines. Since the frequency of 
the 328 cm~^ mode in the Bi compounds does not vary 
with the number of Cu02 layers, we conclude that it 
has 0(2) and 0(3) character. This phonon line is much 
broader than those at 427 cm~^ in Bi2201, and it resem- 
bles the modulation induced splitting, which also sup- 
ports the given assignment. 

A weak structure is present in the T12201 spectrum 
around 230 cm~^, and it could correspond to the 220 
cm~^ mode in Bi compounds. The shape of this mode 
is temperature independent, which is the characteristic 
behavior of in-plane phonons. We therefor attribute this 
mode to an in-plane optical phonon. 

The presence of the 328 cm^-'^ mode with c-axis 
polarization and the absence for this polarization of 
the 220 crp_r|L. mode is in agreement with rigid-ion 
calculationsElH (see the 223 cm^^ mode of T12212 in 



Table [II and the absence of its counterpart in T12201). 
However, these experimental results contraijli£;t the shell- 
model calculations for T12201 (Table |!|),EI where the 
355 cm~^ mode is lacking. 

The low frequency modes are difficult to assign, since 
the heavy atoms involved in this kind of motion have 
comparable masses. The only mode that shifts with the 
change from the single to the double layer is the 114 cm^^ 



mode in Bi2201 (98 cm"! in Bi2212). The tendency to 
shift to a lower frequency is similar to those of the 383 
cm~^ mode, and indicates that the motion of Cu atoms is 
involved. The other two phonons should include the Cu- 
O bending vibrations to minimize interaction with Ca, 
and consequently, motion of Bi, Sr, Tl, and Ba atoms to 
preserve the center of mass. 

In all the compounds there exists a rather large dis- 
crepancy between lattice dynamical calculations and ex- 
periment, as shown in Table |l|: The calculations pre- 
dict the additional mode of 0(2)-0(3) character between 
400 and 500 cm^^ with a rather large oscillator strength. 
Experimentally, this mode seems to be absent. Sec- 
ond, the mode with the largest TO-LO splitting is pre- 
dicted to have lower frequency in T12201 (294-341 cm~^ 
mode) than in T12212 (355-437 cm^^ mode). The con- 
tradiction with the experiment can be readily observed in 
Fig. n. A possible reason for this inconsistency is proba- 
bly the fact that the early lattice dynamical calculations 
have been_adjiisted to data obtained on polycrystalline 
materials, Ej~tfl on which the in-plane and out-of-plane 
optical phonons can not be separated. In fact, in T12223 
and Bi2212 single crystals the in-plane phonons haiie 
been observed at 415 cm"-"^ and 420 cm"-"^, respectively.El 



B. Resonances in intrinsic Josephson junctions 

In the superconducting state the high-Tc compounds 
are considered to form a stack of superconducting Cu02 
planes separated by the non-superconducting heavy el- 
ement rock-salt-type layers, where each pair of Cu02 
planes forms a Josephson junction. The interplane cur- 
rent showing typical Josephson behavior was observed 
in / — y characteri^ics of thin mesoscopic structures of 
Bi2212 and T12223.E3 Subsequent measurements revealed 
a number of small intra-gap-|ifesonances on the resistive 
branch of the I — V curves. e3H The following investiga- 
tions ruled out many known effects related to the super- 
conducting energy gap or geometrical resonances. The 
I — V characteristic of a stack of junctions was shown 
to be a m£re superposition oi I — V curves of individual 
junctionsJ£j which implies that all junctions are indepen- 
dent and that the intra-gap structures occur at character- 
istic bias voltages in every Josephson junction. Recently, 
Helm et alS3 proposed that oscillations of the electric 
field within a junction, excited by the bias voltage due 
to the ac Josephson effect, couple to the IR active lattice 
vibrations and resonate at the longitudinal modes of the 
c-axis phonons. The comparison to the optical data was 
not done, because the resonances lay well within the gap 
and the corresponding frequencies are difhcult to access 
with presently available samples. Moreover, in the usual 
normal incidence technique, the longitudinal frequencies 
are obtained by Kramers-Kronig analysis and can be bi- 
ased by noise and an error in absolute reflectivity at low 
frequencies. 



In Fig. |5| the measured function (1 — i?p)/(l + Rp), 
where maxima correspond to the c-axis longitudinal fre- 
quencies is shown for Bi2212 together with the positions 
of the four pronounced resonances observed in Ref. |5^. 
In the superconducting state the phonon structure does 
not reveal substantial changes as it can be seen from the 
10 K curve. In the lower panel of the Fig. ||, the optical 
data for T12201 and T12212 are compared to the tunnel 
data for T12223. Two phonons in the Bi2212 curve are in 
good agreement with the I — V resonances at 99 and 173 
cm~^, and at least one coincides in the case of T12223. 
The latter phonon is observed in T12201 and T12212, but 
its position is almost unchanged in the systems with one 
and two Cu02 layers per unit cell. We can expect that 
the third layer in T12223 should not substantially shift 
its frequency. 

This agreement let-ius conclude that, as stated in the 
model of Helm et a/.,a the optical phonons play the main 
role in the structure observed in the /— V characteristics. 
However, more resonances are observed in the tunnelling 
experiment. The reason for this discrepancy between the- 
ory and experiment is that Helm et aZ.Ej consider the 
barrier filled with the macroscopic quantity of the mate- 
rial with properties similar to the bulk. In the following 
we will show that the additional resonances can corre- 
spond to zone-boundary modes, which are excited due to 
the non-uniform field distribution in the tunnelling mea- 
surements. Indeed, the electric field E in the interlayer 
Josephson junction oscillates only between two adjacent 
Cu02 planes and vanishes outside the junction. Thus, 
in the wavevector space, -E'((7c) spans the whole Brillouin 
zone following the Fraunhofer diffraction pattern. In op- 
tical measurements, the electromagnetic wave decays in 
the sample over the penetration depth length scale, which 
is much larger than the interplane distance. The corre- 
sponding distribution in qc space is represented by the 6- 
function near qc = Q. The electric field in the Josephson 
junction excites phonons for all qc, and the induced polar- 
ization has peaks at the maxima of the phonon density 
of states. To illustrate this qualitatively we calculated 
the phonon dispersions for the three lowest frequency 
modes using a reduced rigid-ion model. The curves are 
plotted in Fig. ^ja, together with the electric field distri- 
bution E{qc). For this calculation we represented each 
atom in the unit cell by a plane perpendicular to the 
c-axis. The force constants, which include the nearest 
and next-nearest interactions in the unit cell are adopted 
from the Ref. ^. The phonon density of states -D(w) 
is shown in Fig. pp. It can be seen that the polariza- 
tion should peak at the zone-center and zone-boundary 
phonon frequencies. Fig. pd shows the frequency depen- 
dence of the polarization for the uniform field distribu- 
tion (dashed curve) and for the non-uniform box-car dis- 
tribution (solid curve). In the latter case the electric 
field amplitude is finite (and constant) only within a sin- 
gle unit cell as shown in Fig. gc. To simulate this non- 
uniform field distribution, we extended our model to 1000 
unit cells within the chain, with a finite dissipation in- 



troduced to avoid problems due to the circular boundary 
conditions. Fig. |6|d shows that the polarization indeed 
has maxima both at uj{qc — 0) and uj{t:/c). 

The zone-boundary phonons modes can become op- 
tically active if lattice distortions occur in the crystal 
structure. In Fig. g we plotted the data for Bi2201. 
The orthorhombic and monoclinic distortions of 14/ mmm 
structure are stronger in Bi2201 than in Bi2212, and this 
can explain why the 130 cm~^ mode is observed in the 
single layer compound. It should be noted that the Ra- 
man active modes could also be excited in the intrin- 
sic Josephson junction measurements due to the strong 
(~ 10^ V/cm) polarization of the media by the DC bias 
voltage. 



VI. CONCLUSION 

We performed grazing incidence reflectivity mea- 
surements on bismuth and thallium superconducting 
cuprates. The technique was shown to permit accu- 
rate determination of the c-axis longitudinal frequencies. 
We fitted spectra with the Drude-Lorentz model and 
showed that the values of the phonon parameters agree 
with those obtained in normal incidence measurements. 
Therefore, we argue that the grazing incidence technique 
is an excellent tool for probing the c-axis properties in 
high-Tc's on samples with small size in the c-direction, 
which is often a problem in normal incidence measure- 
ments. 

We compared the c-axis TO and LO phonon frequen- 
cies in five structurally similar compounds, Tl2Ba2Cu06, 
Tl2Ba2CaCu208, Bi2Sr2Cu06, Bi2Sr2CaCu208, and 
Bi2_xPbxSr2CaCu208, with lattice dynamical calcula- 
tions and proposed a new assignment for the main oxy- 
gen modes. The major difference is the assignment of 
the 366 cm~^ mode in 2212 compounds, which we con- 
sider to be the in-phase motion of 0(1) and 0(2) against 
Cu. We observed a substantial discrepancy between the 
available lattice dynamical calculations and our experi- 
mental data: The calculations predict a phonon mode 
between 400 and 500 cm^^, which we do not observe. 
Shell model calculations indicate that the oxygen mode 
with the largest TO-LO splitting in the double layer 2212 
compounds should be higher in frequency than those in 
the single layer 2201 compounds, which is again incon- 
sistent with our data. 

The comparison of the LO frequencies obtained by our 
method with those inferred from intrinsic Josephson ef- 
fect measurements shows that the latter technique - due 
to the non-uniform electric field distribution - is sensitive 
both to zone-center and zone-boundary phonons. 
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TABLE I. c-axis phonon parameters of the T12201 and 
T12212 at room temperature. The values between the brack- 
ets were obtained at 8 K, eoo = 4. 

TlaBaaCuOe 



vt ■ cm 



77.9 (77.9) 
156.5 (156.0) 
385.5 (385.7) 
603.2 (603.3) 



vl \cm 



7i 



111.0 (111.3) 

158.0 (157.8) 

427.1 (426.7) 
628.3 (628.5) 



26.5 (26.5) 

8.6 (6.1) 
8.3 (6.5) 

8.7 (8.6) 



S, 



4.5 (4.5) 
0.12 (0.11) 
0.94 (0.95) 
0.28 (0.28) 



TlaBazCaCuaOs 


i>T,- [cm^^] 


i>Li [crn^^] 


7j [cm '] 


Sj 


152.0 (152.0) 

339.0 (339.0) 

598.1 (602.0) 


156.4 (157.0) 
366.3 (367.0) 
614.7 (617.0) 


12.2 (9.0) 
20.8 (18.0) 
17.1(12.0) 


0.13 (0.13) 
1.02 (1.10) 
0.22 (0.23) 



TABLE 11. c-axis phonon parameters of the Bi2201 and 
Bi2212 at 300K (8K), e^o = 4. 



TABLE IIL Calculated values of the c-axis phonon param- 
eters 
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Jia et a/.Nt 
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work 


Kulkarni et alM 
Assignment utj {vl^) 


Assignment 




83.0 


(87.0) 


118.3 


(115.6) 


15.0 (12.0) 


0.35 (0.36) 


!^T, 


144.5 


(148.7) 
(198.3) 


169.7 

217.2 


(175.0) 
(218.0) 


11.0 (11.0) 
16.0 (12.1) 


0.70 (0.68) 
0.68 (0.60) 
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160 
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FIG. 1. Experimental (symbols) and fitted (solid lines) grazing angle reflectivity Rp for the p-polarized light measured at 
80° angle of incidence at the room temperature for (a) bismuth and (b) thallium compounds. Sketched in the inset is the 
schematic diagram of the grazing angle technique. 
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FIG. 2. The c-axis conductivity (a) for the Bi and (b) for Tl compounds, obtained from the fit of the reflectivity to the 
Drude-Lorentz model. 
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FIG. 3. The c-axis loss-function (a) for the Bi and (b) for Tl compounds, obtained from the fit of the reflectivity to the 
Drude-Lorentz model. Shown in (c) is the experimentally measured (1 — R-p)\nab\l(X + Rv) function of T12201, which peaks at 
the longitudinal frequencies. 
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FIG. 4. Crystal structure of Tl2Ba2Cu06 (left) and Tl2Ba2CaCu208 (right). 
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FIG. 6. (a) Simulated phonon dispersion of three branches in the low frequency range for Bi2212 and Qc-dispersion of the 
electric field in tunneling measurements, (b) The phonon density of states corresponding to the dispersion curves in (a), (c) 
Tunnel vs optics electric field distribution within the crystal. In the tunnelling measurements, the electric field is localized 
within one junction, (d) Polarizability for the uniform (dashed line) and non-uniform (solid) electric field distribution, simulated 
on the stack of the afo-planes, whith each plane representing an atom in the unit cell. 
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